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Abstract 

The volumes, spectra and geodesies of a recently constructed infinite family of five- 
dimensional inhomogeneous Einstein metrics on the two bundles over S*^ are exam- 
ined. The metrics are in general of cohomogeneity one but they contain the infinite 
family of homogeneous metrics T^'^. The geodesic flow is shown to be completely in- 
tegrable, in fact both the Hamilton- Jacobi and the Laplace equation separate. As an 
application of these results, we compute the zeta function of the Laplace operator on 
j^pA fQj- large p. We discuss the spectrum of the Lichnerowicz operator on symmet- 
ric transverse tracefree second rank tensor fields, with application to the stability of 
Freund-Rubin compactifications and generalised black holes. 
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1 Introduction 

Compact Riemannian Einstein manifolds may be used as basic building blocks for solu- 
tions to higher dimensional gravity and supergravity theories. An important example in 
recent years has been Preund-Rubin compactification in the context of the AdS/CFT cor- 
respondence. Here one has backgrounds such as AdS^ x supported by 5-form flux. The 
Einstein manifold M encodes geometrically the properties of the dual conformal field theory 
such as the R-symmetry and central charge [1]. Another set of examples are generalised 
D-dimensional black holes, where the horizon is given by an arbitrary Einstein manifold M 
rather than the usual sphere S''^~^[2] . In fact, the two examples we have just given are re- 
lated. A generalised cone spacetime over the Einstein manifold M is a Ricci flat Lorentzian 
metric 

ds^ = -dt^ + dr^ + r'^dslt . (1) 

Generalised black holes may be thought of as black holes in generalised cone spacetimes, see 
(20) below. Instead of black holes, we could have appended three more flat directions and 
considered an extremal D3-brane sitting at the tip of the cone. The near horizon geometry 
of this D3 brane is then AdS^ x [3]. 

Homogeneous Einstein manifolds have been known for some time, well-known examples 
in the physics literature include the round spheres S"^ and the five dimensional T^'* spaces. 
It is harder to find explicit inhomogeneous Einstein metrics on compact manifolds. When a 
metric is found, it is then useful to study its properties, both to achieve a better geometric 
understanding of the manifold and with a view to physical applications. One key question 
for inhomogeneous manifolds is the separability of partial differential equations such as 
the Laplace equation and the Hamilton-Jacobi equation for geodesies. Separability is an 
important first step for being able to perform calculations involving the metrics. 

In four dimensions, the first explicit inhomogeneous compact Einstein metric was con- 
structed by Page [4]. Topologically the manifold is the nontrivial 5^ bundle over 5^, which 
is isomorphic to CP^#CP^. The method of [4] was generalised recently to obtain, amongst 
other results, an infinite series of inhomogeneous Einstein metrics in five dimensions on 
S'^ X and on the nontrivial bundle over [5]. The infinite series is parameterised 
by two integers {k\, k2)- When ki = k2 = k, the metrics become the series of homogeneous 
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metrics T*''^. This construction was then further generahsed to higher dimensions in [6]. In 
this paper we shall be concerned with the properties of the five dimensional metrics labelled 
by ikuk2) [5]. 

There have been two other recent explicit contructions of inhomogeneous Einstein mani- 
folds. Firstly, Bohm has constructed an infinite family of inhomogeneous metrics on S"^ • • • S"^ 
and products of spheres [7]. These metrics can be unwieldy because the metric functions 
are not given explicitly, but as solutions to nonlinear Ordinary Differential Equations. How- 
ever, some properties are known[8]. Secondly, an explicit infinite family of inhomogeneous 
Einstein-Sasaki manifolds were recently constructed on S'^ x in five dimensions [9] and 
also generalised to higher dimensions [10]. 

An important question concerning Preund-Rubin compactifications and generalised black 
hole spacetimes is whether they are stable. In both cases, the question of stability reduces 
to the question of whether the spectrum of the Lichnerowicz operator acting on rank two 
symmetric tensor fields satisfies a certain lower bound [11, 2, 8, 12]. For the AdS/CFT 
correspondence one is more interested in stable spacetimes, as this is when one expects a 
valid duality. For generalised black holes spacetimes, one may also be interested in unstable 
spacetimes because the endpoint of the instability is a nontrivial and interesting dynamical 
question in higher dimensional relativity. Finally, one may also consider the stability of 
generalised black holes with a negative cosmological constant [12]. In this case, the stability 
or instability of the spacetime is related to a poorly understood phase transition in a dual 
thermal field theory induced by the inhomogeneity of the background [12, 13]. 

1.1 Outline of the paper 

The organisation of this paper is as follows. 

In section 2 below, we review the Einstein metrics constructed in [5]. We give a broad 
picture of the discrete moduli space of metrics and we write the metrics in a form where 
the SU{2) X U{1) isometry and the cohomogenity one property is manifest. We then go on 
to characterise the moduli space in terms of volumes and Weyl curvature eigenvalues. The 
latter allows us to comment on the stability and instability of the resulting spacetimes. 

In section 3 we study the Laplacian spectrum on the manifolds. In the homogeneous 
cases, TP'^, we give the spectrum explicitly. The spectrum on T^'i has previously been 
studied in [1]. In the inhomogeneous case we show that the Laplace equation separates 
and for slightly inhomogeneous metrics we give the spectrum as a perturbation about the 
homogeneous cases. 
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In section 4 we use the results on the Laplacian spectrum to calculate the zeta function 
on TP'^ for large p. The zeta function contains information such as the thermodynamics of 
a scalar gas on the Einstein background. Ultimately, one would like to calculate the zeta 
function for an inhomogcncous background, as the corresponding thermodynamics may 
provide some insight into the phase transition expected on backgrounds with a region of 
large curvature [13]. 

In section 5 we study the geodesies of the Einstein metrics. We show that the Hamilton- 
Jacobi equation separates and describe some qualitative features of particle motion. For the 
homogeneous metrics we give a full action-angle analysis of the geodesies. This allows us to 
calculate the periods of geodesies and the semiclassical energy spectrum. The semicalssical 
spectrum agrees with the full spectrum of the Laplacian up to ordering ambiguities in 
quantisation. In the inhomogeneous case the action-angle problem reduces to an evaluation 
of elliptic integrals. These may be calculated in a perturbation about the homogeneous 
metrics, and again the semiclassical energy spectrum has a good agreement with the full 
Laplacian spectrum. 

Section 6 is the conclusion and contains suggestions for future research. 

2 The metrics 

The complete Einstein metrics presented in [5] depend upon two integers ki and k2 which are 
the Chern numbers of a principle bundle over S^. The 5-manifold is then an associated 

bundle over S'^. Topologically there are two such bundles because 7ri(S'0(4)) = Z2. If 
ki + k2 is even the bundle is trivial and if ki + k^ is odd the bundle is non-trivial. 

The metrics may be written as 

dsl = h{efde^ + h{ef {dx^ + sin^ xdr/^) + aij{e){dilj' + cos xdvKdil^ + cos xdv) , (2) 

where the ranges of the angles are < 9 < 7r/2, 0<x<7r, 0<?7<27r and < < 47r/|/sj|. 

The general expressions for the metric quantities h, b, Uij and the integers ki are given 
in appendix A. Throughout we will use a normalisation such that the Ricci scalar of the 
metric is R = 20. The explicit formulae for the metric depend upon two constants i>i and 
f2 which depend in a complicated implicit fashion (given in appendix A) on the integers ki 
and k2- We shall be concerned with the case when ui,U2 > 1. In this case it appears that 
for each pair of positive integers ki, k2 there is a unique pair z^i, 1^2 > 1. Note that there is 
a symmetry interchanging ki and k2 and with it vi and ^2- 
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We will shortly consider the moduli space of Einstien metrics in some detail. Let us first 
exhibit some special cases. A schematic summary of the following statements is contained 
in Figure 1 below. 

• If fci = A;2 = k, hence 1/1 = 1/2 = v, we obtain the homogeneous metrics known in the 
physics literature as T^'^. After the change of variables [5] 

(5 = 26, ^=i(V;2-V^i), i=i(V,i+V;2), (3) 
the metrics takes the standard form for T^'^ 

"^^^"■^ " A{2 + v'^) ^"^^^ ^ ''''' ^"^^'^ ^ 4{2u^ + 1) ^^^^ ^ ''''' ^"^"^'^ 

where k = + 2)/(2i/^ + 1). The case (^1,^2) = (Ij 1) is the Einstein-Sasaki metric 
known as T^'^. 

• The cases (^1,^2) = (0,2) or (2,0) have z/i = 1 or 1^2 = 1 and 1/2 or 1^1 arbitrary, 
respectively. These are on the boundary of the cases we consider. The metric is independent 
of 1/2 or i/i and coincides with the round metric on S^/Z2. 

•. We can let f 2 — 00 with ui finite. One finds that 



We call this the vertical limit. The metric approaches the round metric on S'^/Zfc^. It has 
an orbifold singularity along a circle. This may be described as follows. The metric takes 
the form 

dsl^ = de^ + sin^ 9 dsl + cos^ edip^ , (6) 

where ip has range 2-jt and 

ds| = [(c/^'i + cos xdi])^ + dx^ + sin^ xdv^] ■ (J) 

The angle V'l is identified modulo Therefore, there is an orbifold singularity along the 
circle at ^ = 0. Locally the singularity is M^/Z^^ x S^. 

A word of caution about the vertical limit is necessary. Given that k2 ^ and k2 is an 
integer, in fact the only solution is k2 = 0. Thus, although the limiting orbifold metrics we 
have just described exist, there are no Einstein metrics 'near' the limiting metric. A more 
interesting limit is the limit in which we allow both vi and 1/2 to become large. 
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• For large values of vi and 1^2 we find that 

„,^l-^kf + kl), U2^j-ikl + k^^). (8) 

We call the limit in which both ui and 1^2 become large, with the ratio f 2/1^1 ~ ^1/^2 = Q 
fixed, the rational limit. Near the rational limit, the metric becomes 

J 2 ,n2 cos^ e + q"^ sin^e 3 sin^ ^cos^ e{uji - uj2f + (1 + q^){q'^ sin'' Ouj + cos^ Olo^) 

'^^ + 4(1 + g2) ^^"52+ 4(l + g2)2(cos2 + g2sin2e) 

(9) 

where = dip'^ + cos X'^'?- The actual limiting metric itself is not necessarily Einstein in 
this limit. However, we do have an infinite sequence of Einstein metrics as we approach the 
limit. Therefore, this limit is a richer source of Einstein metrics than the vertical limit. 




Figure 1: Schematic representation of the limits in the moduli space of Einstein metrics. 
The figure is a little misleading in that the rational limit should tend to the upper right 
hand corner, and there are no metrics near the vertical limit. 



2.1 The isometry group 

An SU{2)l X U{1)l x U {1)u isometry of the metric may be made manifest by rewriting the 
metric in terms of left-invariant SU{2) x U(l) forms 

(Ti = cos I Idx + smI jsmxdv, 

(72 = — sin I ax + cos I sin xdrj , 
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•73 = ^ {dipi + dip2) + cos xdr) , 

= \ (#1 - di^2) ■ (10) 

The metric becomes 

dsl = h^de'^ + 6^ [al + aj) + 011(0-3 + a^f + a22{(T3 - cr^f + 2ai2(c73 - (T4)((^3 + 0-4) , (11) 
which we write as 



dsl = h^de^ + {al + al) + {fa^ + a^f + g^aj , (12) 



with 

^2 _ (Q22 - auf ^2 _ 4(a22Qii - 012) ^ _ Qii + 022 + 2ai2 ^^^^ 

On + 022 + 2ai2 ' On + 022 + 2ai2 ' an — 022 

The second form of the metric written here is useful for calculation using the obvious vielbein 

= hd6 , = bai , = 6(72 , = c{f(^3 + C4) , = 5(74 . (14) 

The left-acting isometries follow from writing the metric in terms of left-invariant forms. 
The remaining U{1) symmetry is the usual right-acting ^7(1) isometry of the squashed 
three-sphere. It is present because the metric does not depend on 77. 

2.2 Moduli space of metrics: Volumes and Weyl eigenvalues 

Figure 2 shows the discrete moduli space of Einstein metrics with 1 < z^i,z^2 < 100. Each 
point corresponds to an Einstein metric. 

So far wc have paramctcriscd the moduli space in terms of i^i , U2 or ki,k2- Although 
these quantities were useful for discussing various limits, they do not have an invariant 
geometric meaning. In this subsection we will consider how two geometric properties of the 
metrics vary as we move around the moduli space. The two properties will be the volume 
of the manifold and the eigenvalues of the Weyl tensor. The reason for choosing these 
quantities is that they are related to the physics of generaUsed black holes and Preund- 
Rubin compactifications constructed using these metrics. We will elaborate on physical 
implications in the next subsection. 

The volume of the manifolds is straightforward to calculate and is given by 



VOi^Z^i, Z.2j ^ ^4^, ^ _ ^^2^2^^^2 + ^2^2 _ 2)(^2 + ^2^2 _ 3) " ^'^) 
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Figure 2: Moduli space of Einstein metrics. Crosses correspond to topology x S'^. 
Circles have the topology of the nontrivial bundle over S'^. 

The eigenvalues of the Weyl tensor acting on symmetric tracefree tensors are given by 

Ca^'^hcd = Khab . (16) 

For inhomogeneous manifolds the eigenvalues k = k{x) depend on the position. The max- 
imum value taken by an eigenvalue on the manifold, kq, gives a lower bound for the Lich- 
nerowicz spectrum [8]. In five dimensions and with our normalisation for the Einstein met- 
rics, the bound is 

Al > 4(5 - no) . (17) 

The classical stability of generalised black holes and Freund-Rubin compactifications de- 
pends on the Lichnerowcz spectrum, so kq is an interesting quantity to consider. 

Let us assume that z>2 ^ J^i- Results for the opposite case may be obtained by inter- 
changing Vl and 1^2- In the present case, it turns out that the maximum value of the Weyl 
eigenvalues is achieved at ^ = and may be shown to be given by a fairly simple expression 



(yl - 1)^(2 + 2vl + V9 - 2vl + ^v\) 



(18) 
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To derive this expression, one should work with an orthonormal set of vielbeins (14) and 
consider a basis of fourteen symmetric tracefree matrices. The action of the Weyl tensor 
(16) on this basis will then define a matrix whose eigenvalues are required. 

We can gain some intuition about instabilities by considering the cigcnniodc correspond- 
ing to the eigenvalue (18). At ^ = 0, and using the tangent space basis given by the vielbeins 
(14) we have 

/4 \ 

-3 + x 

-3 + x , (19) 

-2 -2a; 

V W 

where x = 3uj - i/9 - 2i''f + Qvf. One may check that Ca^b^hl^ = Koh^^, aX 9 = 0. We 
expect that for the unstable spacetimes the unstable mode will be concentrated near 6 = 
and will be well approximated by hP at that point. It is curious that the mode does not 
depend on U2- 

Table 1 shows the values of the volume and the maximum Weyl eigenvalue in two limits. 



Limit 


Volume 


Weyl eigenvalue 


Homogeneous: 1/1 = 1/2 = 1^ 


(1 + 1.2)5/223/2^3 


2 + 2i/2 + V9 - 2z/2 + 9i/4 


(2 + Z/2)2(2zy2 + l) 


l + Z/2 


'Rational limit': 

1^2/2^1 = > 1 fixed 


/ 1 TT^ 

V Q 





Table 1: Volumes and maximum of Weyl eigenvalues in various limits. 



We see that in the rational limits, including the homogeneous limit — > 00, the volume 
goes to zero. On the other hand, the Weyl eigenvalue remains finite and may be tuned to 
any value allowed by the integrality constraints: fcj G Z"*". 

2.3 Physical consequences: Entropy and stability 

Five dimensional compact Einstein manifolds appear in a simple way in two contexts in 
higher dimensional gravity. Firstly, they can be used to construct generalised black hole 
spacetimes in seven dimensions 

ds"" = -f{r)de + ^ + r'^dsl (20) 
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with /(r) = 1 — (1 + r^/L^)(r+/r)^ + r^/L^. This is a solution to Einstein's equations, with 
a possible negative cosmological constant —l/L?. The horizon is at r = r^. The volume 
of the five dimensional manifold becomes the area of the event horizon and is therefore 
proportional to the entropy of the black holes. 

The classical stability of generalised black holes has been studied recently, both with 
vanishing cosmological constant [2, 8] and with a negative cosmological constant [12]. The 
latter case has interesting field theory implications using the AdS/CFT correspondence [13]. 

For generalised black holes there is a simple criterion for classical stability. The criterion 
depends on the minimum Lichnerowicz eigenvalue of the horizon. For a five dimensional 
horizon, and with vanishing cosmological constant, the criterion is [2] 

Al > 4 <^ stability. (21) 

Therefore, from (17) we see that if any of the metrics have kq < 4, they will result in stable 
spacetimes. 

For large generalised black holes with a negative cosmological constant, r+/L >> 1, the 
stability criterion is [12] 

> -A^ X ^ stability, (22) 

where is a positive 0(1) coefficient that may be determined numerically [12]. In order 
to be unstable, the maximum Weyl eigenvalue kq must therefore be very large. 

A second type of solution that uses Einstein manifolds are Freund-Rubin compactifica- 
tions with a 5 form field strength 

1/2 

vol(M^) . (23) 

Remarkably the stability of these spacetimes [11, 8] is according to precisely the same 
criterion as for generalised black holes (21). This may be due to the relation between 

generalised black holes and Frcund- Rubin conipactifications pointed out in the introduction. 

The relation between the Lichnerowicz spectrum and the eigenvalues of the Weyl tensor 
(17) is such that it allows us to prove the stability of the above spacetimes, but not insta- 
bility. However, it was found in [8] that large positive values of kq tend to suggest unstable 
spacetimes. 

One can check from (18) that kq > 4 for all values except z/i = z/2 = = 1, which has 
precisely kq = 4. Therefore the only rigorous statement about stability we can make is that 



8(D - 2) 
D-6 
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the u = 1 metric gives marginally stable spacetimes. This metric is non other than T ' so 

we have reproduced a known fact [14]. However, it seems very likely that all the remaining 
metrics will give unstable spacetimes. The instability means that these metrics have limited 
interest for the AdS/CFT correspondence, but give rise to interesting physics of generalised 
black holes. 

Like the Bohm metrics [8, 7], the moduli space here includes regions where the maximum 
Weyl eigenvalue is becoming arbitrarily large. For example, q is arbitrary in Table 1. 
This suggests that the lower bound on the Lichnerowicz spectrum may become arbitrarily 
negative. This in turn implies that the metrics give not only unstable flat space black holes, 
but also unstable large black holes in anti-de Sitter space according to the criterion (22). 

Unlike the Bohm metrics however, the metrics we consider here are given explicitly, 
without need of numerical calculations. Therefore it is substantially easier to study quantum 
fields on these backgrounds. We begin this study below by computing the zeta function of 
a free scalar field in some cases. Ultimately, one would like to understand better the nature 
of the field theory instability that seems to appear when the maximum Weyl eigenvalue 
becomes large [13]. 



3 Laplacian spectrum 



In this section we show how the Laplace eigenvalue equation for the metrics we are consid- 
ering separates. In the homogeneous cases we find the spectrum explicitly. The equation 
we need to solve is 



with 



where 



and 



□ = + 



d . d 1 . ,2 ^ ^ 



b{9y sin X ^ ^dx'^ b{6)'^ 
1 ( d 



63 



sinx \dr] 
1 



cosx 



d 



d d 
+ 



h{e)b{e)'^^deta{e) oe 



6(^)Vdeta(^) d 



80 



where a*-' is the inverse matrix of Uij. 



(24) 
(25) 
(26) 
(27) 



The eigenvalue equation can be separated by using the left acting SU{2) x U{1) isometry 



11 



to rewrite the Laplacian in terms of the symmetry generators 



= 

6 = 
6 = 

^4 = 



— cot X sin 

— cot X cos 

d d 

+ 



2 



d d 

+ 



■ d d 
+ 



+ cos 



— sm 



■0^ + -0^ 
2 



^ sm 

sin X 



cos (V^) 5 



smx 



d d 



Define the SU (2) quadratic Casimir 



(28) 



sinxdx 

The Laplacian may thus be written as 



d . d 
smx^i h 



d d 
+ 



-\ 2 



dx dtjj^ dtp'^ 



(29) 
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Now note that ^^,.^3,^4 all commute with each other. Wc know their eigenvalues from 
SU{2) group theory. Wc must be careful to account for the coordinate ranges correctly. 
Let ^ be a simultaneous eigenfunction. Then 



^[6-^4] 



= —inicf) , 
= -in2<j), 

-i(i + i)0, 



n2 e ^z, 

J G Vu{0} 



(31) 



The range of ni+n2 is restricted by —j < ni+n2 < j ■ All the eigenvalues have a degeneracy 
of 2[7] + 1, where [x] denotes the integer part of x. The allowed values of ni + 712 follow 
from the standard angular momentum result ni + 712 G {—j, — j + 1, — j + 2, . . . , j} and then 
retaining only the values for ni + n2 that are consistent with (31). We state the allowed 
values more explicitly below in (42) for the homogeneous case. 

One is left with an ordinary differential equation for the 6 dependence of the eigenstate 



jjj + 1) - ("1 + "-2) 



-(f) + a^-'niUjcf) = X(f) . 



(32) 



Solutions of this equation will introduce another discrete parameter I labelling the eigen- 
values. Write the eigenvalues of the full Laplacian as Xj^ni,n2,l- The zeta function for the 
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Laplacian on the Einstein manifolds is thus 



c-(«) = EEE 



2b1 + i 

A? 



(33) 



with the summations restricted by (31) and the comments following. 



3.1 Homogeneous metrics: spectrum 

In the homogeneous case one has i/i = U2 = v and 

ki = k2 = k 



(34) 



2z/2 + 1 ■ 

The Laplacian simplifies. The relevant terms are collected in appendix A. The equation 
becomes 



1 



d 



sin 9 cos 6 dO 



sin 9 cos 9 



d_ 
89 



+ 4 



{2u^ + 1) 



(Z^2 + 2)2 



z/2 sin2 ( 



+ 



Z^2 cos2 I 



-4 



(2^.2 + i; 



j-(i + i)'^ + 6 



(2Z.2 + 1] 



(m + n2f(t> + 



(l + r^^ 



(35) 



1.2 + 2 ■'^■^ ' ' "(1/2 + 2)2^'"^ ' '-^^ ' 1.2 + 2 
This equation may be solved in terms of hypergeometric functions. If one rewrites the 
equation as 



d 



sin 9 cos 9 89 



sin V cos ( 



89 



+ 



^2 



+ 



sin2 cos2 9 



A0, 



(36) 



then two linearly independent solutions are 

(j)^ = cos^^^sin^ 92F] 



, /l + A±S- (l + A)i/2 1 + A±S + (1 + A)V2 
^ 2 ' 2 



± S; cos^ 6* 



(37) 

Regularity at = 7r/2 requires that ±5 = Regularity as — requires both ^ = 
1^1 > and further that the hypergeometric function be polynomial. We will see below that 
the solutions are in fact Jacobi polynomials. This condition for (37) to be a polynomial is 
that 

l + |y4| + |B|±(l + A)V2 



-L. 



lg Z+U{0}. 



The solution of this equation is A = + 1) with 

2Z.2 + 1 



l = L + 



viv^ + 2] 



{\ni\ + |n2|) 



The eigenvalues of the Laplacian are thus seen to be 



(38) 



(39) 



(40) 



A 



'j,ni,n2,( 



Zy2 + 1 



(i.2 + l)(z.2 + 2) 



(ni + 712)^ . 



In the following subsection we clarify the values that n^, j, I may take. We will recover the 
spectrum (40) from a semiclassical quantisation of geodesic energies in section 5.2 below. 
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3.2 Towards the zeta function 

To calculate the zeta function later it will be convenient to parameterise the eigenvalues 
with the following four integers 

J,LgZ+U{0}, iVi,iV2GZ. (41) 

The eigenvalues (40) are then calculated using 

m = -Ni , 

j = ^\Ni+N2\ + J, 

I = ^(|iVi| + |iV2|) + L. (42) 
The zeta function may be written 

<-oW-EE E t (43) 

L=0 J=0 Ari=-oo Ar2=-oo Ni,N2,J,L 

In the following section we will use this parameterisation to calculate explicitly the zeta 
function for the homogeneous metrics as u ^ oo. These are the metrics T^'^ with p large. 

3.3 Inhomogeneous metrics 

The differential equation for the 6 dependent part of the solutions (32) is significantly more 
complicated in the general inhomogeneous case. The first step towards extracting informa- 
tion from the equation is to clarify the structure of the singular points of the equation. 
Let us write z = cos 29. Then (32) becomes 

^0 + P(z)^</. + g(z)<^ = u;(z)A0, (44) 

where 

P{z) = + + , (45) 



and 



r^t \ t \\ ^1 I ^2 , Qs .R1.R2.R3 
Q(z)-a;(z)A=^^-^ + ^^^ + ^^-^ + ^ + ^ + ^, (46) 



with zq = 2 cos ^0 — 1 and 
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The pole zq arose as a root of = in [5] . The coefficients Qi are given by 
and 2 

^3 ^2-uf-ui){i~ufui)y h + k2 J ■ ^^^^ 

The other coefficients Ri satisfy a relation Ri + R2 + R3 = since Q{z) = o(l/z^). The 
explicit form of these coefficients is somewhat complicated, so we leave them to appendix B. 
Thus (44) is a Fuchs type differential equation with 4 regular singular points including 00. 
Such equations have a canonical form studied by Heun. We transform (44) to the canonical 
form in appendix E. 

Prom (45) and (48), the eigenfunction around z = 1 {0 = 0) takes the form 

<P{z) = {l-zt^/''^\f{z) (50) 

with f{z) an analytic function at z = 1, while the eigenfunction around z = —1 {9 = 7r/2) 
takes the form 

4){z) = {l + z)\''^'^^\g{z) (51) 

with g{z) an analytic function at z = —1. In the homogeneous case, we saw that the regular 
function is given by a polynomial. 

At infinity, the regular singular point implies that solution behaves as 



0(2) = z-^±V^-Qi-Q^-Q3+i^o-i)Ri+{^o+i)R2^^^ (52) 

where h{z) is analytic at infinity. For the homogeneous metrics — 00, and (52) needs 
to be recalculated. However, infinity remains a regular singular point in the homogeneous 
limit. 

3.4 Slightly inhomogeneous metrics: spectrum 

In this subsection, we compute the Laplacian spectrum for metrics that are slightly inhomo- 
geneous. We can do this using Quantum Mechanical perturbation theory. In the rational 
limit at least, there arc metrics sufficiently near the homogeneous metrics in space for 
which perturbation theory is applicable. 
The perturbation we consider is 

— = 1 + e , with e < 1 . (53) 
^1 



15 



Equation (44) may be multiplied by — 1 so that it is expressed in terms of a self-adjoint 
operator H. We will check self-adjointness shortly. 



i7<^ = - 1) + P{z)j-^m + Q{z)^{z) 



= Cj{z)\<t>{z) , (54) 
with uj{z) = (z^ — \)u}{z). The coefficients are then perturbed to first order as 

P = p(0) + ,p(i) ^ g = Q(0) + gg(i) ^ = + ed;^ . (55) 

The zeroth order operator is thus 

= {z^ - 1) + + ■ (56) 

Using the explicit form of the function 



the operator becomes 
which is manifestly self-adjoint with respect to the inner product 



Hi^) = ±y-l)± + ^,^-l)Qio)^ (58) 



<<f>l,(f)2>= J dz(f)l{z)(f)2{z). 



(59) 



Expanding the eigenfunction (f> = (f>^'^^ + e<l)^^^ and eigenvalue A = A'^^^ -|- eA^^^ we have the 
well known result from first order perturbation theory 

,(1) ^ <p^"\(//W-A(")^n-)),,(0)> 

<</.(o),w(o)0(o) > ■ ^ ' 

In the present case we have 

ij(l)_A(0)^W = (^2_i^|^p(l)^+g(l)_;^(0)^(l)^ 

^ .^(1 + .^) , d 

(z/2-l)(zy2+2)^ 

+ Q?'^^+Q^^^'^^ + iR?+4')^ + R?-li'\ (61) 

the actual expressions for R^^^ and Q^^^ are still rather large, so we relegate them to appendix 
B. The denominator in (60) is given by 

^^'^ = 7^1^^ • (62) 



4(2 + z/2) 
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In calculating (61) one should keep ni, n2, j, I fixed under the perturbation 1^2/1^1 = 1+e. 
This is because these arc properties of the zeroth order solution, not of the full equation. 
The relationship between ni,j,l and Ni,J,L is kept to be that of the homogeneous case 
(42) and is not modified. 

To evaluate (60), note that (p^^\z) is just the homogeneous solution that we found 
previously (37). It is convenient at this point to rewrite these solutions in terms of Jacobi 
polynomials. It turns out that, up to an overall normalisation 

0(o)(^) = (1 _ ^)i^ii/2(i + ,^m/2pm\,m)^^^ ^^3^ 

We have used the integers introduced in (41). We recall that the Jacobi polynomials have 
the following definition 

^i^'^H^) = ^(l-^)-"(l + ^)-^(^)"[(l-^)"+"(l + ^)^+'*] 

oc F{-n,l + n + a + /3,l + a;{l- z)/2) . (64) 

The advantage of introducing Jacobi polynomials is that they satisfy identities that will 
enable us to perform the various integrals needed in the evaluation of (60). We see that 
we will need to evalute < 0^, Z(/>(o) >, < l/(z - 1)0(°) >, < (p^°\l/{z + l)(/)(°) > and 
< (/.(°),(22_ 1)^/^^0(0) AH of these expressions may be computed, the required identities 
involving Jacobi polynomials are given in appendix C. 
The result is 

<<^(o)^(^(i)_A(o)cdW)0(o) > 



1.2(1 + 1.2 



(z/2-l)(r/2 + 2) 



m + m + 1 + 2L ^(i) liVil + |iV2| + 1 + 2L ^(i) ^ ^g^^ 



with the coefficients 

^1 



|7Vi| ^1 |7V2 



(2L + |7Vi| + |iV2|)(2L + |Ari| + |7V2| + 2) ' 



2L+|Ari| + |7V2| ^ ' 

This expression does not appear to simplify to a pleasant expression in terms of n^, j, I. We 
will see in section 5.4 below how a very similar result for the perturbed spectrum may be 
found from a semiclassical quantisation of the classical geodesic energies. 
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4 Zeta function for T^'^ at large p 



We will consider the homogeneous metrics in the limit — oo, corresponding to A; — oo 
also. The advantage of taking this limit is that summations of discrete eigenvalues may be 
approximated as integrals over continuous parameters. This fact will enable us to compute 
the zeta function. There are infinitely many metrics in this limit, evenly spaced in v. See 
Figure 2 above. Towards the end of this section we will consider one physical application 
of the zeta function. 

The spectrum becomes (40) 



A = 8j{j + 1) + 4Z(Z + 1) - l|(ni + n2f . 



(67) 



One should express the spectrum in term of the integers of (41) 



rii = -Ni , 
3 = J+^lATi + ATsI, 

I = L + ^{\N,\ + \N2\). 

where J, L G Z+ U {0} and iVj G Z. The degeneracy of the eigenvalues is 2[j] + 1. 
In the large v limit the spectrum may further be simplified to 



(68) 



A = 8 



J+-\Ni+N2\ 



+ 4 



^ + ^(l^l| + l^2|) + l 



(69) 



This step does not hold when A^i + = 0, we will consider this case separately later. Until 
specified otherwise, assume that Ni + N2 ^ 0. The degeneracy may now be taken to be 
2J+||A^i+iV2|. 

The zeta function at large u may therefore be written 

y {2l + l\N, + N2\) 

NML (2 [2i + i|iVi + N2\]' + ^[L + i(|iVi| + |iV2|)] [L + i(|iVi| + [N^D + l])' 

(70) 

Note that is large and J only appears as J/iy. We can approximate this sum by an integral 
when Ni + N2 ^ 0. That is to say, we write 



j=o 



V^^ 



(71) 



The integral is straightforward to perform and the result is 



,2-2s 



8(1 -s) ^ V2 

^ ^ Ni,N2,L 



L + l(|iVi| + |iV2|) 



^+ 2(1^1! + 1^2!) + 1 



l-s 



(72) 
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In order to do these sums, the modulus signs force us to consider four cases separately. 
First note that we can take A'^i > and the remaining eigenvalues have an extra degeneracy 
of 2, which we will include at the end. The four cases then depend on the value of A''2. The 
cases arc: (I) N2 > 0. (II) N2 = -Ni. (Ill) -00 < N2 < -Ni. (IV) -Ni < N2 < 0. The 
total zeta function will be the sum of these contributions 



Cis) = 2 [Ci{s) + Cii{s) + Onis) + Civis)] . (73) 



A check of our result is given in appendix D. 

4.1 Case I: N2>0 

In this case we may write M = Ni + N2 > 0. We consider the case M = later. The sum 
is expressed entirely in terms of M. This introduces a degeneracy of M + 1: the number of 
ways of writing M as the sum of two positive integers. The zeta function becomes 

We may now convert the L summation into an integral for the same reasons as before. One 
obtains 

This integral can be done using, for example, the Maple program. The answer may be 
expressed in terms of hyper geometric functions. Taking the large u limit, one obtains 

£ ^!!^,3-2.£(lziZ^M3-2^(l + M) . (76) 

M=l ^ ' 

But now the remaining sum is just a Riemann zeta function, so the final answer for this 
case is 

Ois) = ^-^^^^-'-''^^^^^f^ K«(2^ - 3) + C«(2. - 4)] . (77) 

4.2 Case II: N2 = -A^i 

In doing this case we will include the A^^i = iV2 = case. In these cases, the integral method 
is not valid. The zeta function contribution from these values is 
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As usual, the summation should exclude the zero mode. This sum may be simplified by 
considering Q = L + Ni. It may be possible to do this sum using contour integration. This 
will not be necessary here; the only point that is relevant is that to leading order there is 
no u dependence and that the first correction is order This term will not give the 

dominant contribution to physical quantities which will be at e.g. s = 0, —1/2. 

4.3 Case III: -oo < N2 < -Ni 

Define the quantities M = Ni + N2 < Q and N = Ni - N2 > 0. The zeta function 
contribution for this case is 

" ■^ + -^(L + iV/2)(L + iV/2 + l) 



^ ^ N,M,L 



(79) 



What is the range of the summation for M and A^? If wc fix M, then we see that N takes 
the values = |M|, |M| + 2, |M| +4, .... Therefore L + N/2 takes the values \M\/2, \M\/2 + 
1, \M\/2 + 2,... with degeneracy L + N/2 + 1 - \M\/2. Defining x = {L + N/2)/l' we may 
rewrite the summation over L and A'' as an integral over x: 

. ^ x+l/u-\M\/2u , 

CiiAs) = — — > / ■ — jdx . (80) 

8(^-l)MtllW(Mi+4x(x + l/.))^" 

One may again do this integral. It is similar to the one considered previously. The final 
result for the zeta function is 

-2^ v^-''%r{2s - 4) [{s - 1)2^1 (2, s,b-s- -1) + {s- 4)2^1(1, g - 1, 4 - g; -1)] 
128 (s-l)(s-2)(s-3)(s-4) 

(81) 

Note that this expression always has one power more in p than the contribution from case I 
(77). Therefore the contribution from case I is negligible unless the expression (81) vanishes. 
But if this expression were to vanish then we would need to calculate the subleading term. 
This will not be a problem in practice. 

4.4 Case IV: -N^ < N2 < 

Define the quantities M = Ni + N2 > Q and N = Ni - N2 > 0. The zeta function 
contribution for this case is 

,2-2s 



M2 4 ^ 
— + -(L + N/2){L + N/2 + 1) 



(82) 
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This is the same expression as before (79) except that now if we fix M the range of N is 
M + 2, M + 4, .... Now L + N/2 takes the values M/2 + 1, M/2 + 2, with degeneracy 
L + N/2 — M/2. Defining x = (L + N/2)/i' as before the zeta function becomes 

To leading order in v this integral turns out to be the same as in the previous case 

C,w{s) = C,nAs). (84) 

4.5 Physical applications of the zeta function 

The most immediate application of the zeta function is to calculate the thermodynamics 
of a free scalar field on the compact manifold M. We can add a mass or coupling to the 
Ricci scalar without changing the zeta function that we have calculated to leading order in 
z/, as can be seen from (69). We denote this coupling with k. More precisely, k will only be 
relevant when |A^i + A'^2| = 0, and we saw that this case gives subleading contributions to 
the zeta function. 

The free energy is given by the logarithm of the partition function on M x S""*^, where 
the has length /? and the temperature of the scalar radiation is T = 1//3 as usual, 

where we have introduced an arbitrary mass scale /x so that the dimensionalities are correct. 
Zeta function regularisation [15, 16] then gives a finite expression for the formal determinant 

^ = -^Ca(0), (86) 

where we have set A = [—{d/dt)'^ — □ + k]///^. The zeta function is analytic at the origin, 
so the free energy is finite. 

We can calculate the free energy in the limits of low and high temperature. In both 
cases, it is convenient to use the following expression for the zeta function 



= ^Indet ' ' \^ ^ , (85) 



1 



00 



U{s) = TiA-' = I dTT^-^Tve--^ . (87) 
i \S) Jo 



The eigenvalues of (d/dt)"^ on the factor are 



27tV 



^n= (yj ri\ neZ. (88) 
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It follows that we may write the zeta function as 

1 °° /"OO ( '2-n\^ 2 

Us) = wT E / dTT'-\-'\^^) " TV^e--^^ , (89) 

n=-oo -^0 

where As refers to the operator acting on the d dimensional spatial section M. 

First consider the theory at low temperature. In this regime we may approximate the 
sum over the modes by an integral 

CO / „ \ 2 

(2 ^ 

2^7_oo"""' - (47r)V2 rV2 " 



^ O-rr I A'n-\l 2 ^1/2 ^ 



n=— oo 



The zeta function may thus be written as 

_ Pfi r(.-i/2) 

- (47r)V2 r(.) ^^-(^-V2)- (91) 

The technique for calculating the high temperature behaviour was introduced by [17]. 
For the terms in the sum with n 7^ 0, it turns out that we can use the Schwinger-de Witt 
expansion for the heat kernel Tr^e^'^'^'^ as r — >^ 0. This will work when the temperature is 
larger than the curvatures, because in this case the dominant contribution to the integral 
comes from small values of r. The Schwinger-de Witt asymptotic expansion is 



00 ^ 

TV.e-^^« ~ V ^r*-'^/^ = -J^ 



<=f/2 + _| 



aoT , 

/X2 



(92) 



where the first few a2k are given in appendix D. If we let denote a generic curvature 
scalar with mass dimension 2k, then we have that a2k = C(£ft*). 

Substituting the expansion (92) into the expression for the zeta function (89) and doing 
the integral over r, one obtains 

r) d °° / /I 2 \ d/2—s—k 

= ^^^^'^ + (4^0^ ?o "^^^ ^ ' " ^ " ^''^ 

where Cr is the Riemann zeta function. We see that the series expansion will be valid if 
I 0'2kP'^^ 1*^ li that is, if the temperature is large compared with curvature scales | $H |^ T^. 

Thermodynamic quantities may be calculated using standard formulae such as = 
d(/3F)/df3 and S = —dF/dT. The results in the low and high temperature limits are 
shown in Table 2. In this table we restore the dependence on the scalar curvature R for 
completeness. At low temperature, the leading contribution comes from C,in{s) and C,iv{s)- 
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Quantity 


Low Temperature 


High Temperature 


f 


H 

2419200 


Mbv 21601/ 


E 


VlOiii/S 
2419200 ^ 


2^/27^6 (6k + i?) V27r4 . 
189z/ 720i/ 


S 


+ --- 


4V27r6^5 (6/^ + i^)^/27r4^3 
3151/ 540z/ 



Table 2: Thermodynamics of a scalar field on M x 5 at high and low temperature. 



5 Hamilton-Jacobi equation and geodesies 
5.1 Separability of the Hamilton-Jacobi equation 

In this section we show that the Hamihon-Jacobi equation for the geodesies of the metrics 
we are considering is integrable. Hence we obtain first order equations for the geodesies. 
The Lagrangian for a free particle is 

C = gabdx^'dx^ = h^e'^ + b^{x^ + sin^ x??^) + aiji'^' + cos xv){ip^ + cos xv) , (94) 

where a dot denotes differentiation with respect to some time parametcrisation of the 
geodesies. The geodesic equations may be separated by using a Hamilton-Jacobi method. 
In this description, the dual momenta are 

The coordinates '0* and rj are cyclic, giving three conserved momenta, p^i = Ji and Prj = J3. 
The first order equations for these coordinates are easily seen to be 

. _ 1 ( Ji + J2) cos X - J3 

2 62 gjj^2 ^ ' 

1/^* = ^a'' Jj — cos xv ■ (96) 
To find the remaining geodesic equations, consider the Hamilton-Jacobi equation 

with E a constant. This equation may be separated by taking 

S = W{e) + U{x) + J^ri + Ji^i + JsV'^ • (98) 
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One then obtains the equations 



+ —^{h-cosx[Ji + J2]f = L', (99) 



and 

1 (dwV ^ ^ L2 

-I- nP.hJ, + — 



h?\de ) ' ' 62 

where L is a constant. Now use the relations 



+ a'^JiJj + ^ = E, (100) 



dU 


dS 


dC 


dx 


dx 


dx 


dW 


dS 


_ dC 




'ae 


~ ae 



2h^e, (101) 

to obtain equations for x ^'Hd 9 from (99) and (100) respectively. These equations can then 
be used to eliminate the time dependence and hence obtain equations for the orbits. For 
example 

h? fdeV - L'^ /b^ - a'UiJi . , , 



62 \dx) L2 sin2 x-\h- cos x( Ji + J^y? 
Similar equations may be derived for -0* and r) in terms of %. In fact, one obtains a simple 
equation for dx/dr] 

This equation describes the projection of the geodesic onto the (x, if) two-sphcrc. One can 
recognise (103) as describing the motion of a charged particle on a sphere with a magnetic 
monopole background. This is not surprising given that the full metric was constructed 
from a principle bundle over the S"^ . It is well known that charged particle motion on 
a sphere with a magnetic monopole background results in closed circular orbits. A little 
algebra shows that the motion (103) describes a circle on the two sphere about an axis n 
with opening angle G. One finds 

and 

n 



j2 J 

1 — -|- sin^ [sin S e^; + cos S e^] + — sin 06^ , (105) 
where S is an arbitrary angle and {ea;,ej,,e2} is the standard Cartesian basis. 
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5.2 Action-angle analysis of the homogeneous case: semiclassical spec- 
trum 

In this subsection we carry out a full action-angle variable analysis of the geodesies of the 
homogeneous metrics. This analysis allows us to compute the frequencies of closed geodesies 
and also allows a semiclassical quantisation of the system. We see that the semicalssical 
spectrum essentially agrees with the spectrum of the Laplacian that we calculated previ- 
ously. 

We work in phase space 

M = T*{S'^ X S^), (106) 

with coordinates (x") = (/? = 29,x,^,'(pi,tp2) and their dual momenta The canonical 
one form a and corresponding symplectic form u are as always 

a=PaAdx'^, uj = da. (107) 

The Hamiltonian for the system is 

H{x,p)=g^''paPb. (108) 

We may write H explicitly as 

^ = ^ 1^+ J2 ^ + fc2 +P^if - 2(PV2 +Pi^i){PH -i?^i) cos/3 + (p^^ - Pi>if)^ 

+ 1 + 1,2 [Px + ~ ^(^^2 +P^i)PvCOSX + (PV2 +P^l) ) 1 

^^^^^"'^ (P^2+P^J^ (109) 



(l + Z/2)(2 + i.2) 

This Hamiltonian system is integrable since there exist 5 = dimM/2 independent functions 
Fa (a = 1 ... 5) such that 

• {Fa,H}^ = 0, 

• {Fa,Fb}o, = 0. 

The functions may be taken to be 

?2 2 1 



-^1 = P(3 + psin2^ ((PV'2 +PVi) - '^{Pi'2+Pi>i)iPi^2-P^i ) cos l3 + (p^^ - p^^) ) 

^2 = Px + ~~ir- (Pv ~ "^(P^i + PV-i )Pv COS X + (PV2 + P^l f) 
^ sm X 

F3 = Pr,, Fi=p^^, F5=p^^. (110) 
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In which case the Hamiltonian (109) may be written 



In fact, the quadratic conserved quantities Fi and F2 are related to two reducible Staeckel- 
KiUing tensors, iiT?'', of the background 

Fi = Kfp,p, (i = l,2). (112) 

This follows from writing the quantities as 

where we have used the generators of the SU (2) x SU (2) symmetry, 

^1 = -cotxsin(('0i +V'2)/2)(pvi +cos((V'i + V'2)/2)Px + ^!^i^^!tlljMI^p^ ^ 

Sin ^ 

^2 = - cot xcos((V'i + 'i/'2)/2)(pv'i - sin(('i/'i + '^2)/'2)Px + !^!!ii^!l±i!^lZ^p^ ^ 

and 

= - cot /3 sin(A;(V'i + V'2)/2)(pv'i + ^^2)/^ + cos(A;('0i + ^2)/2)p/3 
sin(A;(V'i + V'2)/2) 



-V'2)/2)(pv-i+: 
-(pvi -PV'2)A> 



sin/3 

= - cot acos(A;(V'i + V'2)/2)(pvi +^^2)/^ " sin(A;(V'i + V'2)/2)p/3 
cos(fe(V'i + V'2)/2) 



sin/3 

6 = {Pi,i +P^2 )/k , (115) 

which satisfy the relations = -eijfe6 and = -eijkik- 

The five constants of motion allow us to consider the level set 

Mf = {{x,p) : Fa{x,p) = fa}. (116) 

General theorems show that Mf is diffeomorphic to the 5 dimensional torus. The action 
variables are constructed by considering 5 independent cycles in the torus, Ca, and writing 

Ia = ^f a. (117) 



2vr ./c„ 



The definition is in fact independent of the cycle chosen to represent a given homology 
class because the symplectic form cj = da vanishes when restricted to My. From their 
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definition, the variables will be invertible functions of the constants {fa} only. If we 
invert this relationship, we may consider fail)- The key step in action-angle analysis is 
then to note that the Hamilton-Jacobi function may be considered as a function of x"' and 
/„: S = S{x,I). Thus one may define new coordinates 

*° = ^. (118) 

These are the angle variables. The action-angle variables satisfy three important properties. 
Firstly, (0", /„) are canonical coordinates on the phase space 

ui = dcp" A dia . (119) 

Secondly, the variables (f)"' are indeed angles on the cycles Ca 

d4P = 2-K51 . (120) 



ICa 

Thirdly, the time evolution equations for {(j)"-,Ia) are trivial 

ia = 0, 

= dH{I) ^ ^ = j^«mi + const. (121) 

The next step is to choose 5 cycles in Mf. Three cases are particularly straightfor- 
ward. Consider cycles that have tangent vectors and The action variables are 
respectively, using (110), 



h = f3, h = \fi (i = 4,5). 



(122) 



The remaining two integrals are slightly more complicated. They are most easily computed 
by re-expressing the action variable as an integral over a surface in phase space with Ca = 

dSa 

27Tla = j> a= / da= / dpbAdx''. (123) 

First we describe the curves. The two cycles will be taken to be the intersection of Mf with 
the surface generated by tangent vectors (^,^) and (^,^) respectively. From (110) 
we can calculate the equation for the cycles, in both cases the curve takes the general form 

C = {{x,y) e [0,7r] xR\y^ = a^ -{b^ -2bccosx + c^)/sm^x}. (124) 

If the following two conditions hold 

- \bc\ > 0, (a^ - b'^){a'^ - c^) > , (125) 
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then the equation for the curve becomes 

= a^(cosx — Si){S2 — cosx)/ sin^ x, — 1 < (5i < ^2 < 1 , 



with 



1 



^1,2 = ^{bc± V(a2-62)(a2-c2)) 



(126) 



(127) 



The curve C is closed and the area A enclosed by C is 

= |a|7r (2 - V(l + 5i)(l + ,52) - V(l - 5i)(l - <52)) 

= 7r(2|a| - |6 + c| - |5-c|). 



(128) 



By (123) we see that the action variables will be given by the area inside the curve in the 
{x,p) plane. Now consider the two cases: 

(A) Ii{f): In the equation for the curve (124) we have a = f\,b = {f4 + /s)/^, c = 
if 5 ~ f4:)/k- We can take /i > without loss of generality. Then, |-fi(/)| = Ai/27r, where 
Ai is the area given by (128) 



l^i|=/i- 



(129) 



(B) hif)'- Here we have a = f2, b = + f^, c = /a, and taken /2 > without loss of 
generality. Then, |/2(/)| = ^2/27r, where A2 is the area given by (128): 



I/2I = /2 - ^(1/4 + /5 + /3I + 1/4 + /5 - /al) . 



(130) 



Here, we have assumed the condition (125) for /„, which is equivalent to that of the existence 
of periodic orbits in M. 

We can now calculate the frequencies of the orbits. The frequency Jl" is (121) 

dHil) 



din. 



where 



1 + 1/2 

Using (129)and (130) we have 



1 + 1/2 



2(l + z/2)(l + 2z/2)' 



(131) 



(132) 



/i = |/i| + 2(1^41 + I/5I), 

/2 = |/2| + ^(|fc(/4 + /5) + 2/3| + |fe(/4 + /5)-2/3| 



(133) 
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This is a classical version of the Laplace spectrum (40) written with the action coordinates. 

We may use the expression for the Hamiltonian in terms of the action variables (132) 
to perform a semiclassical quantisation of the spectrum. The semiclassical prescription is 
to replace action variables by integers. Consider the following quantisation 

h^L, h^M, h^Ni, h^N2, 
h^J+^\Ni + N2\ - ^{\k{Nr + N2) + 2M| + |A;(iVi + N2) - 2M|) . (134) 

The slightly awkward quantisation of I2 is in order to make contact with our previous 
notation. It is not difficult to see that the above expression sets I2 to be an integer and J 
to be a positive integer. 

Using the definitions of (42) we see that the semiclassical spectrum may be written as 

4(2^^ + 1) 2 , 4(1.^ + 2) 2 + 1) I . ns^l 

This is precisely the same as the spectrum for the Laplacian that we calculated previously in 
(40), up to the usual semiclassical error — j(_7 + 1), — Z(Z+ 1). This agreement suggests 
an approach for calculating the Laplacian spectrum in the inhomogeneous cases. We can 
understand the semiclassical disagreement by comparing a naive quantisation of the classical 
Hamiltonian (109), — > idxo.-, with the Laplacian (25). The quantised Hamiltonian has 
terms 5| and 5^ whilst the Laplacian has 1/ sin/3(?^[sin/35/3] and 1/ sin [sin x^^^] . This 
represents an ordering ambiguity in quantisation and explains versus j(j + 1). 

Finally, we can write down an explicit solution for the functions M^(/3, /) and U{x,I) 
that appeared in the separation of the Hamilton- Jacobi equation (98). Using the following 
formula, which is the integral of (128) in an indefinite form, 

- V(l+«(l.«arc.a.(^/|±||L^) 
, ^,1 - area. (^1*^) 



2arctan W^— (-1 < (5i < ^2 < 1), 



we find that 




Wi(3, 1) = /i(/)K(cos /3; 61 (/), (52(1)), 

U{x,I) = f2{I)K{cosxMl)Ml)), (136) 
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where the parameters 6i{I) and Sj{I) are determined as a function of la by (127) with 
the values for a, b, c given in (A) and (B) above, respectively. W and U are multivalued 
functions which have multiplicities NAi, NA2, with N £ Z. 

5.3 Action-angle analysis of the inhomogeneous case 



After setting z = cos 20, hence pg = —l-pz^X — z'^, the Hamiltonian is given by 

^ ^ ^^h'^{z) ^z) (^^ si^^^^^ ~ ^^^"^^ +P^i)Pv^^^X + iPi>2 + Pi>if cos^ X)^ 

+ a'^{z)p^,p^. . (137) 

This is also an intcgrablc Hamiltonian system. The following functions are mutually com- 
muting conserved quantities, 

Fi = H, 

^2 = P^ + ^^(.Pl-'^{P■>P2+P^l)Pv^O^X + {P^P2+P^plf): 

^ sm X 

F3 = Pt], F4=Pi,i, F5=PiP2- (138) 
As in the homogeneous case, the new momentum coordinates are introduced by 

^(/) = ^ / 

Four of the action variables are the same as in the homogeneous case 

I/2I = /2 - ^(1/4 + /5 + /3I + 1/4 + /5 - /3I) , (140) 

and 

/3 = /3, h = ^f4, h = ^h. (141) 
The coordinate Ii is harder to calculate. On the level set 

Mf = {{x,p):Faix,p) = fa}, (142) 

the z-component Pz is written as 



where 



Q{z) = 7 TT9 + 7 TTo + 7 7 7 ) (144) 

{Z-I)^ {z+lf {Z-Zof Z-1 Z + 1 Z-ZQ ^ ' 
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and the coefHcients Qi,Rj are given by equations (48) and (49) and the equations in ap- 
pendix B, together with the replacements 

+ ni,n2^f4,h, A^/i. (145) 

Given the expression (143) for pz, we would like to calculate the action variable 

h = ^jpzdz. (146) 

Let us introduce 

^{z) = {z- \)\z + lf{z - zofQiz) . (147) 
The function ip is a polynomial of degree 4 with leading coefficient — /i/4 and 

ip{l) = AQi{zo - if < 0, <^(-l) = 4Q2(zo + 1)' < 0, ¥^(^0) = Q3(^o - 1)' < 0. (148) 

The turning points of the geodesies are given by the roots of ip{z). We will consider only the 
case where the polynomial ip has 4 distinct real roots, say (i = 1, 2, 3, 4). Indeed, for large 
values of vi and 1^2, we can find real roots such that either (a) — 1 < cki < 02 < 1 < 03 < 04, 
or (6) cci < 0:2 < — 1 < 0:3 < 04 < 1, depending on the values of /„. Thus the expression for 
the action variable becomes the following elliptic integral 



■^Jai [Z - 1){Z + 1}{Z - Zo) 



(149) 



where ctj and aj are the roots between -1 and +1. 

It seems that real roots ai satisfying (a) or (b) exist in the general case, although we 
have not proved this. One can check that real roots ai satisfying (a) or (b) exist for small 
values of /4, /s and any z/j. Note that roots of the form — 1 < ai < a2 < eta < 1 < 04 or 
ai < — 1 < 02 < CK3 < CK4 < 1 are forbidden by (148). 

As a more explicit example, consider the rational limit z^i, f2 00 with 1^2/1^1 = q > 1 
fixed. Keeping /i, /2, /4, fixed in the limit one finds Qi = 0, (1 = 1, 2, 3) and 

R, = {l + q^)ll-^l, i?2 = -(l + g2)^ + |., R^ = ^R^-R^^ (150) 

This limit is given by (192) and (193) in appendix B with TVj = 0. Thus the polynomial 
becomes 

viz) = iz-l)\z + lfiz-zofQ{z) 

^ -^{z-l){z + l){z-zo){z-zi), (151) 
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where zq = {q^ + l)/{q'^ — 1) and zi = zo{l — 8/|//i). Positivity of (p and \z\ < 1 requires 
one of four cases 

(/) -1< z < zi<l< zo, {II) -1 < z < 1 < zi < Zo, 
(III) zo<zi<-l<z<l, {IV) zo < -K zi < z < 1. (152) 

Near the rational Hmit, the 4 real roots {— 1, 1, zq, -zi} change as 

-1^-1 + di, 1^1 + 52, zo^z^, zi^zl, (153) 

where Zq,zI are the new values of zq, zi and Si are evaluated up to the order u^^ as 

8Q2(l + ^o) ^ 8Qi(l-zo) 

fl{l+Zl) fl{l-Zl) 

This implies the inequality {a) — 1 < ai < a2 < 1 < < a4 or {b) ai < a2 < —1 < < 
a4 < 1 by Qi < and /i > 0: 

(/),(//)^(a), {III),{IV)^{b). (155) 

Therefore the action variable is of the form (149). The integral may be evaluated in the 
four cases 

{1) -1 < z < zi < l< Zq : 



with 



|7i| = ^^\ (n(7r/2, (^1 + l)/2,,.) - F{7r/2,K)) , (156) 

■iry/2{zo - zi) 



^=J^-^^^^^. (157) 
V 2(zo-^i) ^ ' 



(II) -l< Z <1< Zi< Zq: 



\h\ = .f'^'' =n(7r/2, 2/(zi + 1), «) , (158) 

^^^/{zo - l){zi + 1) 



with 



2{zo - zi) 

{zo-i){zi + iy 



(159) 



(III) Zo < zi < -1< z < 1 : 



III I = ^^(^1 + 1) n(7r/2,2/(l-zi),K), (160) 

7rV(^o + i)(-2;i - 1) 



with 



. = J , '^']~''\, . (161) 
V (^o + l)(^i-l) ^ 
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(IV) zq < -l< zi < z <1 : 

, , , ^(zi + 1) 



ny/2{zi - Zo] 

with 



(n(7r/2, (1 - zi)/2, k) - F{7t/2, «)) , (162) 



.= X-^^^P^. (163) 

In these expressions F and 11 denote complete ehiptic integrals of the first and third kind 
respectively. 

5.4 Slightly inhomogeneous metrics: semiclassical spectrum 

In this subsection we calculate the action variable for a small perturbation about the homo- 
geneous metrics. This is the classical computation corresponding to the spectral calculation 
of section 3.4 above. 

The homogeneous metrics are given by oo. In this limit 

ip{z) ^ zi 0{z) + o(zo) , (164) 

which represents the transition from the Riemann surface = ^p{z) of genus 1 to = 0{z) 
of genus 0. We therefore have 

where is a polynomial of degree 2: 

^{z) = Qi{z + if + Q2{z - 1)2 + (i?i - R2){z - l){z + 1). (166) 

The coefficients Qi and Ri are given by (48), (49) and (188), again with the replacements 
of (145) 

J(J + 1)^/1, ni,n2^h,h, A^/i. (167) 
The action variable is calculated using the integral (128) to give 



|/i(/)| = vTO + IQ^^^RTTR^-vW- VTO- (168) 

Solving this expression for fi = H recovers the energy given in (132). 
Now consider a perturbation away from homogeneity as in section 3.4 

— = 1 + e , with e < 1 . (169) 
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In this case we have a polynomial of degree 3 by ignoring the terms that are o(e^), 

0{z) = Qi{z + 1)2 + Q2{z- 1)2 + - R2){z -l){z + l) + z{z -l)iz+ l){Ri + R2). (170) 

The coefficients Qi and Ri are given by (48), (49) and (183), (184) with 1^2/1^1 = 1 + e. The 
last term in the polynomial (170) has coefficient Ri + R2 which is order e, since R1 + R2 = 
for e = 0. The integral for the action variable can be evaluated as a perturbation of the 
homogeneous limit. We find 

|/i(/)| = ^\Qi\ + \Q2\-Ri + R2-VW\-V\Q^\ 



2(|/i| + ^Q^+ VW)' 

Expressing the coefficients in terms of the homogeneous quantities /4,/5, namely Qi = 
-11/4: + eQf^ and Q2 = -^H^ + eQ^^\ together with (189) and (190) in appendix B, we 
have 

homogeneous i -j- Z/ 



h = P. 



where 



= q'? ^Qf -Rf -m.\^\H^\h\){^^'^ 



\H \h\ 



(ii - iim? + R?) 



(173) 



(2|/i| + |74| + |/5|)2 • 

The perturbed energy is just H = fi. Therefore the semiclassical spectrum is computed 
by setting — > L, 74 — Ni and I5 — N2 in equations (172) and (173). It is satisfying 
to see that the resulting semiclassical spectrum agrees with the Laplacian spectrum (65) 
within a semi-classical approximation. This provides a check on our calculations. We can 
understand why the first term in (65) is lacking from the semiclassical result by again 
comparing a naive quantisation of the classical Hamiltonian (137), Pa id/dx"", with the 
Laplacian (25). There is an ordering ambiguity, and we see that the P{z)d/dz term in the 
Fuchs equation does not exist in the quantiation of the classical Hamiltonian. 



6 Conclusions and future directions 

We have studied the Laplacian spectrum, Lichncrowicz spectrum and geodesies on an infinite 
family of five dimensional inhomogeneous Einstein metrics. The moduli space of metrics 
included an infinite sequence of homogeneous metrics, T^'^. For the homogeneous metrics, 
we were able to give very explicit results for the Laplacian spectrum and for the frequences 
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of closed geodesies. Further, we were able to use the explicit spectrum on T^' to calculate 
the zeta function on these metrics at large p. 

The inhomogeneous metrics are harder to study. We have shown that the Laplace 
equation and Hamilton-Jacobi equation may be separated for these cases. Wc found some 
pcrturbative results for the Laplacian spectrum for slightly inhomogeneous spaces. However, 
it seems that the full inhomogeneous spectrum will require numerical calculations or more 
sophisticated methods than we have used. 

For the geodesies in the inhomogeneous case we have identified four of the five action 
variables. The calculation of the remaining action variable reduces to an elliptic integral 
with an underlying elliptic curve of genus one. If one could calculate the remaining action 
variable, then it is possible that the semiclassical quantisation of the system would give 
some insight into the Laplacian spectrum for the inhomogeneous metrics. Semiclassical 
quantisation of the homogeneous and slightly inhomogeneous cases gives a good agreement 
with the Laplacian spectrum. 

The Lichnerowicz spectrum contains information about the stability of Freund-Rubin 
compactifications and generalised black holes constructed from Einstein metrics. We saw 
that of the family of metrics we considered, only T^'^ can be shown to give a stable Freund- 
Rubin compactification and generalised black holes with vanishing cosmological constant. 
We suspect that the remaining spacetimes are all unstable. Unstable generalised black holes 
are interesting as it is unclear what the endpoint of the instability will be and the instability 
may give rise to interesting dynamics, analogous to what has been discovered recently for 
the black string instability [18]. 

We saw that the moduli space contains metrics with an arbitrarily large maximum of 
Wcyl eigenvalues. This occurs in the rational limit at large inhomogencity q, sec Table 
1 above. This fact suggests that the minimum Lichnerowicz eigenvalue may become ar- 
bitrarily negative [8], giving rise to unstable generalised Anti-de Sitter black holes. The 
instability of these black holes predicts a thermal instability of a dual theory propagating 
on the corresponding Einstein metric [13, 12]. If one could calculate the zeta function for 
these backgrounds, the corresponding thermodynamics may help elucidate the nature of the 
predicted field theory instability. Perhaps the Laplacian spectrum in the large-g rational 
limit is calculable? One should note that it is not certain that the dual instability will 
necessarily be present for a free scalar field, as duality relates the AdS-j black hole to the 
the strongly coupled thermal theory living on M5 branes which is certainly much more 
complicated. However, it seems possible that a thermodynamic instability for field theories 
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on a curved background with regions of large curvature is a generic phenomenon. 

Finally, it would be interesting to perform an analysis similar to ours for other Einstein 
metrics. In particular, Einstein-Sasaki metrics always give stable Freund-Rubin compactifi- 
cations and stable generalised black holes [8] . Therefore a study of the Laplacian spectrum 
of the five dimensional Einstein-Sasaki metrics constructed in [9] may have interesting ap- 
plications to the AdS/CFT correspondence. These metrics are also cohomogeneity one and 
it seems clear that the equations can also be separated in these cases. 
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A Metric terms 



We normalise the metrics so that the Ricci scalar is i? = 20. 



h\9 
aii{9 
a22{0 
ai2{e 



with 



(1 



cos 



z/n sm 



//f cos^ 6 — sin^ ' 



- ulf{2 - vl - - ulul) (1 - nl cos2 - i/| sin^ 9) sin^ 9 



4(1 



.If {2 



J^l - - ^1^2) (1-2^1 COS 



- vf cos^ 9 — sin^ 9 

2 ...2 a _ ^^2 g-j^2 pQg2 



4(1 + ^,4^2 + ^2^4 _ 3^2^2)2 



v{ cos^ 



sin 



sm 



' cos 



4(1 + ufvl + ulu^ 



1 — l^f COS^ I 



— Uo sm 



(1 - ufu^){l - uf cos^ 6* - z/| sin^ 9) 



,2, A 



4(1 + lyfiyi + vlvi 



2 _ vl{\ - vlvl) 



2 _ vl{\ - vlvl) 



The metric has the following integrality conditions 

- ^2) (2 - A - 



ki 



ko = 



1 + ufu^ + vfu'2 — 3z^^z^2 
U2{l-uf){2-u!-u!ul) 



G Z. 



1 -I- ufiy^ + - 2,vlvl 



2,, 4 



,2,, 2 



(174) 
(175) 
(176) 
(177) 
(178) 

(179) 

(180) 
(181) 
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In the homogeneous case, 1/1=1/2 = 1^, 



a 



(z/2 + 2)(l + l/2) V I/2C0S2^ 
12 _ 2(2i/2 + 1)2 

(l/2 + 2)(l + Z/2) ' 

^ li^^^jLil^ h~'^ = ^ ^l , Vdetaa sin6'cos6'. (182) 
1 + i/^ 1 + i/^ 

B Coefficients in the inhomogeneous equation 

The coefficients i?j are given by 

+ 1) - (ni + n2f)B Air+nl , ^2(1 - z/i^)ni , A^jl - uf)nin2 
-n-l — :r7:: ^ :^ 1 tttz r 



8(1 -/X?) 2 8(1-/^2) 4(1-/^2)2 

(l-^D(i_^2^2)^ 



8(2-i.2_^,2)(l_^2) 



2\ ' 



(183) 



D (i(i + 1) - (ni + n2)2)B ^2r-ni Ai(l - i/|)nf ^12(1 - z/|)nin2 
^2 = tt:; H 



+ 



8(1-^2) 2 8(1 -/xi) 4(1 -/.2)2 

(l-z.|)(l-i.H)A 



8(2-I.2_^2)(i_^2)' 

R3 = -R1-R2, (184) 



where 



A2 
A12 



4(1 + 


ufu^ + i/f - 


31.2^1)' 


Z.2(l 


-^2)2(2-^? 


- ^1)^ 


4(1 + 


uful + l/2j,4 _ 


3^M)2 


^1(1 


-^.2)2(2- ^2 


- ^1)^ 


4(1 + 


ufu^ + i^f 1^1 - 


3^.^1)2 




(2 - ^2 _ ^2) 


3 


4(1 + 


1/^1/2 + Z/2jy4 _ 





and 



Cl{l-vlf , C2(l-^.?) C3(l-i^?) (1-^.2)2 



fi ..2\3 0/1 ..2\5 



(l-/.2)3 2(1-^2)2 2(1 -/.2)2 4(1 -/.2)2' 

ci{i-4f C2{i-ul) c^{\-vl) {\-yl? 

(l-/.2)3 2(1 -/.2)2 2(l-;.2)2 4(l-;.2)2' 



(186) 
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with 



ci 

C3 



- - ujul) 
2(1 - vl - vl) ' 

2 

-2v\ + 1^2+1^1 + I'li'l - vlv'. 



2(2 - ul 
2vl - vl - i/| - ulul + vful 
= 2(2 -.f-. I) • ^'''^ 



In the homogeneous case, ui = 1^2 = we have 



Ri = —R2 



il+2u^){j{j + l)-{ni + n2f) {l + 2u^f{nl + nl) (1 + 2r/2) ^712 (l + r^^)A 
2(2 + 1/2) + 4zy2(2 + zy2) + 2(2 + 1/2)2 8(2 + 1/2)' 

i?3 = 0. (188) 



We also collect here the terms for the perturbation about the homogeneous metrics. 
The required terms for Q^^^ are 

(1) _ 4(2z/2 + l)(t/4 + 2z/2 + 3) 2 

^1 " (i.2 + 2)3(z/2_l) '^1' 

0« - 2(2z/2 + l)(7^^ + 3^2^2) 

^2 (^2 + 2)3(z/2-l)^2 "2, (189) 

whilst the terms for R^^^ are 

i/^(-l + 3^2^^4) . , 20 + 56^/2 + 37r/^ + 8^/6 - 4^/8 

^1 - (^2_l)(^2 + 2)2 (j(j + l)-(^l+"2) ) 4(^2 _l)(^2 + 2)3 

4 +101/2 + 1 lz/^ + 16i/6 + 4z/^ 2 i^2(-6 - 5i/2 + lOz/"^ + 4z/6) 

^ 4l/2(z/2-l)(z/2 + 2)2 "2+ 2(z/2-l)(z/2 + 2)3 ""^""^ 

+ ^! A(o) 

^ 4(z/2-l)(l/2 + 2)2 ' 

^2 - (^2_l)(^2 + 2)2^^^-^ + ^^ (m + n2))+ 4(^2 _l)(^2 + 2)2 

8 + 24i/2 + 32z/^ + 41i/6 + 16i/® - 4z/^° 2 z^2(6 + 7z/2 - 8z/^ + 41/^) 

4z/2(i/2-l)(z.2 + 2)3 2(z/2-l)(i/2+2)3 ""^""^ 

+ -:rA(°). (190) 

4(l/2 - l)(l/2 + 2)2 ^ 

Finally, the coefficients simplify in the rational limit, i/i,z/2 — >■ oo with 1/2/2/1 = g > 1 
fixed. One obtains 

.0 = ^ . (191) 
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The coefficients (46) become, introducing the integers Ni by Ui = kiNi/2, 

^^ = -^^^- ^^ = -^^^- Q3 = -^(iVi+A^2)^. (192) 
Furthermore, 



2 



^2 = -(1+.^)^+^ 



^ + i(2 + - - q')Ni + T^(l + 2g2)iV,iV2, 



i?3 = -R^-R^. (193) 

C Properties of Jacobi polynomials 

All the following identities involving Jacobi polynomials may be found, for example, in [19] 
section 8.96. 

The first equation we use allows us to express zPji''^\z) in terms of pji^'^\z), P^^\z) 
andp(^f(z) 

2(n + l)(n + a + /3 + l)(2n + a + /?)P^+f (z) 

= {2n + a + /3+l)[{2n + a + p){2n + a + p + 2)z + a^ -p^]Pj^"''^\z) 

-2{n + a){n + P){2n + a + p + 2)P^"f (z) . (194) 

The following equation allows us to express (1 — z'^)dP^'^\z) / dz in terms of P^'^\z), 
zPt^\z) ^ud pI^^J^Xz) 

i2n + a + p){l-z^)-^Pt^\z) 

= n[{a-P)-{2n + a + /3)z]Pi"'^) (z) +2{n + a){n + /3)Pi"'f (z) . (195) 

Once we have used the above relations, we use two integration results. The first is 

(1 - ^)"(1 + zfpj^'l^\z)Pi"'f^\z)dz 

= , if m 7^ n 

2«+/3+ir(Q, + 1)T{0 + n + 1) 



n!(a + P + 1 + 2n)r{a + p + n + 1) ' 
The second relation we will use is 



if m = n . (196) 



7-1 n\ a L [a + p + n + 1) 
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D Check of the zeta function 



The zeta function on an n dimensional compact manifold (in our case n = 5) , 

C(^) = E . (198) 

i;Ai>0 

can be shown to converge absolutely in the region *Kes > n/2, and can be analytically 
extended to a meromorphic function of s in the whole complex plane. The poles are located 
at s = n/2 — k, {k = 0, 1, 2, • • •) and their residues are given by 

^^^-/^-^= (4.)n/2r(n/2-fc) - ^'''^ 
If the operator is A = — □ + k then the first few are given by 

ao = vol(M), 02 = -ao + \ i dxR , (200) 

and 

Is our computation of the homogeneous zeta function able to reproduce any of these 
results? We have just recalled that there should be a pole at s = 5/2 with residue 

vol(M) , , 

Res.-5/2 = 7 — / , , . (202) 

(47r)V2r(5/2) ^ ^ 

In the large u limit, this gives 

a/2 

Res,=5/2 = ^ . (203) 

We can compare this with our results. We have already noted that case I is always subleading 
compared to cases III and IV. We see that cases III and IV do have a pole at s = 5/2. The 
total residue of ({s) in equation (73) turns out to be 

Res,=5/2 = ^ ■ (204) 

Which is precisely as required! We see that case II cannot contribute to this expression 
because it has no terms that are 0{l/i'). Thus we have a rather nontrivial check on our 
calculation. 

It is not possible to check the other residues because we have only worked to leading 
order in v. The other poles do not appear to this order. 
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E Heun's equation 



In order to transform the Fuchsian equation of section 3.3, equations (44) to (46), to a 
canonical form write 

</.(z) = (1 - z)V^\l + z)V^^{z - zo)^^/(-2) . (205) 

The new variable f{x) with x = {1 — z)/2 satisfies Heun's equation 



where 



dx^^ '''' \x ^ X — 1 ^ X — xq J dx^ ~^ x{x — l)(x — xq) ' ^ ^ 



a + /9 = 2 + 1^ V\Q7\^ , 



a(3 = 2\Y,,/\Q~\+Y^^J\QiQj\ + {Ri + R2)xo-R2\ , 

\i=l i<j 

7 = l + 2vlQ^, 



S = l + 2vTO> 
e = 1 + 2^0^, 

q = {V\Q^\ + V\Q^\+'^V\QiQ2\)xo + yAQ^\ + 2y/\QiQ3\ + 2Rixd207) 



One has 



l + a + /3-7-(5-e = 0. (208) 



It is known that Huen's equation admits an expression in terms of elliptic functions and 

this expression is closely related to the Inozemtsev system. For example sec [21, 20]. From 
these references, if all of 7, 6, e are half-odd-integer one can obtain exact solutions of Huen's 
equation. Unfortunately in our case this condition is not satisfied: 7, 5 are integers and 
e G K. 
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